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We study analytically the steady-state regime of a net- 
work of n error-prone self-replicating templates forming an 
asymmetric hypercycle and its error tail. We show that the 
existence of a master template with a higher non-catalyzed 
self-replicative productivity, a, than the error tail ensures the 
stability of chains in which m < n — 1 templates coexist with 
the master species. The stability of these chains against the 
error tail is guaranteed for catalytic coupling strengths (K) 
of order of a. We find that the hypercycle becomes more sta- 
ble than the chains only for K of order of a 2 . Furthermore, 
we show that the minimal replication accuracy per template 
needed to maintain the hy percycle, the so-called error thresh- 
old, vanishes like \J n/K for large K and n < 4. 

PACS: 87.10.+e, 87.23.Kg 



I. INTRODUCTION 

The limitation of the length of a genome by the repli- 
cation accuracy per nucleotide (q) has led to a deadlock 
in the theories of the origin of life based on the evolution 
of competing self-replicating polynucleotides. According 
to Eigen's quasispecies model which may serve as 

a paradigm here, polynucleotides have to replicate with 
high accuracy in order to reach a certain length, a re- 
quirement that is impossible to fulfill without the aid 
of specialized catalysts. However, to build those cata- 
lysts a blueprint is necessary that amounts to a large 
genome (the nucleotide sequence), which itself cannot 
be maintained without the catalysts. In particular, for 
polynucleotides of fixed length L, the quasispecies model 
predicts the existence of a minimal replication accuracy 
per genome Q c = q^, below which the genetic informa- 
tion is irreversibly lost. This information crisis has been 
termed error threshold transition. Above Q c the popula- 
tion is composed of a master copy together with a cloud of 
structurally similar mutants (quasispecies) |l],|2| . Equally 
important is the finding that, except in a trivially degen- 
erate case, two or more quasispecies cannot coexist 
thus precluding the coexistence of templates (i.e. polynu- 
cleotides) sufficiently different from each other to code for 
any useful set of catalysts. Although it has been claimed 
that the information crisis is not really a fundamental 
issue, since the error threshold transition appears only in 



some pathological, discontinuous replication landscapes 
Q , the coexistence problem seems to be more pervasive, 
as it is associated to the form of the growth functions in 
the chemical kinetics equations [^J. 

In order to circumvent the aforementioned limitations 
of the quasispecies model, Eigen and Schuster proposed 
the hypercycle 0, that is, a catalytic feedback network 
whereby each template helps in the replication of the next 
one, in a regulatory cycle closing on itself. This model 
has gain plausibility when the ability of polynucleotides 
to help propagate each other was established experimen- 
tally through the study of the catalytic activity of the 
RNA (ribozymes) |8|,[|. Interestingly, though the error 
threshold phenomenon has traditionally been considered 
the main motivation for the proposal of the hypercycle 
(see 10 1, for instance), most of the seminal works in this 
field have dealt with the coexistence issue only, as they 
assume perfect replication accuracy for the hypercycle 
elements |7],[ll| . In this case an arbitrary number of tem- 
plates permanently coexist in a dynamical equilibrium 
state; if n > 4, however, the template concentrations 
vary with time j?| , periodically decreasing to very small 
values. In practice, large hypercycles are therefore sus- 
ceptible to extinction via fluctuations, see e.g. hence 
the information gain due to the coexistence of different 
templates in the hypercycle may not be very impressive 
after all. Furthermore, we will argue in this paper that 
coexistence in the absence of a stable equilibrium can 
also be achieved by a simpler arrangement, namely, the 
free chains, in which the cyclic order of the catalysts is 
interrupted. 

The effect of error-prone replication (mutation) in the 
hypercyclic organization was investigated by introducing 
a mutation field as a perturbation of the error-free ki- 
netic equations This approach, however, is not very 
appropriate to study the error threshold phenomenon, 
since the results obtained cannot be easily compared with 
those of the quasispecies model. In this sense, a better 
approach is to assume the existence of a special class of 
templates with no catalytic activity, so-called error-tail, 
that appear as a consequence of the replication errors of 
the hypercycle elements p~^-|Tq| . However, the particu- 
lar catalytic network investigated extensively within that 
framework was not the hypercycle, except for a short dis- 
cussion in fTi| , but the fully connected network in which 
each element helps the replication of all other elements 
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of the network |15|,^6[. (Clearly, in the case of n = 2 ele- 
ments, these two networks become identical Such a 
network is more robust than the hypercycle since the mal- 
functioning or extinction of one of its elements does not 
compromise the whole network. Nevertheless, besides its 
aesthetic appeal, the cyclic coupling of the hypercycle 
seems to be more realistic 0. 

The goal of this paper is to investigate analytically the 
steady-states of a deterministic system comprised of two 
parts, namely, a hypercycle made up of n self-replicating 
templates I±, I2, ■ ■ ■ , I n an d its error tail I e . These parts 
are coupled such that any erroneous copy of the hyper- 
cycle elements will belong to the error tail. The focus of 
the present analysis is on the location in the parameters 
space of the model (i.e. replication accuracy per tem- 
plate, non-catalyzed and catalyzed productivity values, 
and hypercycle size) of the regions of stability of the di- 
verse possibilities of coexistence between the templates 
composing the hypercycle. In particular, we give empha- 
sis to the characterization of the critical parameters at 
which the hypercycle becomes unstable against the error 
tail. 

The remainder of the paper is organized as follows. In 
Sec. [n] we present the chemical kinetics equations that 
govern the time evolution of the system and motivate 
the specific choice of the parameters used throughout the 
paper. The fixed-points of th e kinetic equations are ob- 
tained anal ytically in Sec. Ill and their stability discussed 
in Sec. IV, The phase-diagrams showing the regions of 



stability of the diverse coexistence states are presented 
and analyzed in Sec. |v|. Finally, some concluding re- 



marks are presented in Sec. VI 




FIG. 1. The system composed of a hypercycle of size 
n = 3 and its error tail 7 e . The thin arrows represent 
the non-catalyzed self-replication reactions and the thick ar- 
rows represent the self-replication catalytically assisted by the 
neighbor template 



where xq = x n and 

n 

* = X i ( A i + K i x i-l) + A e 



(3) 



1=1 



is a dilution flux that keeps the total concentration con- 
stant, i.e., ii + x e — 0. As usual, the dot denotes a 



II. THE MODEL 



time derivative. Henceforth we will assume that 15 1 



We consider a system composed of a hypercycle made 
up of n elements I\ , . . . , J„ and its error tail I e , as il- 
lustrated in Fig. [I]. In contrast to the so-called ele- 
mentary hypercycle we assume that the templates 
are capable of self-replication with productivity values 
Ai (i = l,...,n) and A e . Moreover, as usual, the growth 
promotion of template Ii as a result of the catalysis from 
template 7j_j is measured by the kinetic constants Ki. 
The key ingredient in the modeling is that in both pro- 
cesses of growth of template Ii the probability of success 
is given by the parameter Q £ [0,1], so that an erro- 
neous copy, which will then belong to the error tail, is 
produced with probability 1 — Q. Hence the concentra- 
tions Xi (i = 1, . . . , n) of the hypercycle elements and the 
concentration x e of the error-tail evolve in time according 
to the kinetic equations 



±i = Xi (AiQ + KiXi-xQ - $) i=l,...,n (1) 



and 



(Ae-ty + il-Q^XiiAi+KiXi^) (2) 



y^xj+x e = i. 



(4) 



i=l 



Clearly, this formulation is equivalent to considering 
polynucleotides of length L — > 00 whose replication accu- 
racy per nucleotide q goes to 1 such that the replication 
accuracy per genome is finite, i.e. q L — > Q. In this 
limit, the back-mutations from the error-tail elements to 
the templates that compose the hypercycle, as well as 
the mutations between those templates, can be safely ne- 
glected. Hence, mutations can only increase the concen- 
tration of the elements in the error-tail. The advantage 
of working in this limit is that the error threshold transi- 
tion can be precisely located by determining the value of 
Q at which the concentration of a relevant template van- 
ishes. For finite L, as well as for finite population sizes, 
the characterization of this transition is more involved, 
being achieved through the use of finite-size scaling tech- 
niques [ pj| . 

In this work we consider the single-sharp-peak replica- 
tion landscape |l|,Q, in which we ascribe the productiv- 
ity value Ai = a > 1 to the so-called master template 
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I\ and Aj = A e = 1 to the n — 1 other elements of the 
hypercycle as well as to the error-tail. Also, for the sake 
of simplicity we set Ki = K for all i. The motivation 
for this particular choice of parameters is the observa- 
tion that the emergence of the hypercycle requires both 
spatial and temporal coexistence of the templates form- 
ing the network, and this can be achieved by a quasis- 
pecies distribution, which guarantees the coexistence of 
the master template and its close mutants, despite the 
purely competitive character of the quasispecies model 
p9[ . Once the coexistence is established, the appear- 
ance of catalytic couplings between the templates is not 
a very unlike event. Of course, as soon as those cooper- 
ative couplings become sufficiently strong to balance the 
competition imposed by the constant concentration con- 
straint, the mutants will certainly depart from the master 
template due to the relentless pressure of the mutations, 
so that no trace will remain of the original quasispecies 
distribution. 



III. FIXED POINTS 

Let us distinguish between surviving templates Xi > 
and extinct templates Xi = 0. A survivor Ij is said 
isolated if Xj-\ = Xj+\ = 0. Hence, 



ij = xj (Q - $) j > 1 



(5) 



and 



i B = s e (l-*) + (l-Q) 



l) + l+K £ ^-i 



(6) 



In the steady-state regime, Eq. (||) yields $ = Q which, 
for Q < 1, is incompatible with Eq. (|^), since the term 
within brackets in this equation is positive. Therefore, all 
isolated survivors with the exception of the master tem- 
plate are unstable against the error tail. Next consider 
the following chain of surviving templates: 



Xi{Q - $) 



x i+ i = x i+ i(Q + KQxt - $) 



Xi+2 = x i+2 {Q + KQx i+ i - $) 



(7) 



(8) 



(9) 



0, i.e., Xi — 0. So, there is no fixed point corresponding 
to such a chain. Any chain of survivors therefore must 
start with template n or 1. In the first case we get $ = Q 
from x n = and aQ + KQx n = Q from i\ — 0, yielding 
x n = (1 — a)/K < 0, which rules out this possibility. 
The equilibria of interest for our study thus are cither 
the interior equilibrium in which all templates survive, 
or a fixed point that corresponds to a chain of survivors 
beginning with I\ . 

Accordingly, we define a m-coexistence state as the n- 
component template vector x = (x±, x%, . . . ,x n ) in which 
the first m components are strictly positive and the rest 
are equal to zero. Clearly, given the template vector x 
the concentration of the error tail x e is determined by the 
constraint (^). In the following we solve analytically the 
kinetic equations in the steady-state regime ii = Vi. 

The simplest fixed point is the zero-coexistence state 
(to = 0) which corresponds to the solution x\ = ... = 
x n = and x e = 1, existing for the complete range of 
parameter values. 

In the case of chains, i.e. < to < n, the steady-state 
solutions of Eqns. ([!]) and (|2|) are straightforward. In 
fact, since x n = by definition, we get $ = aQ from 
±i = which then yields 



a- 1 



x\ = x 2 = ■ ■ ■ = x m -i = 



K 



Next we insert this result in Eq. (||) to obtain 

Qa-l . a — 1 

(m — IJ 



a - 1 



K 



(12) 



(13) 



However, since Xi £ (0, 1) Vi, this solution is physically 
meaningful in the region K > a — 1 and Q > Q m where 



1 1 / \ / >9 

= - + — (m-l)(a-l) 2 . 
a Ka 



(14) 



We note that the 1-coexistence state (quasispecies) is 
obtained by setting m = 1 in Eq. ( |l3| ) and its region 
of existence is simply Q > 1/a, since the other condi- 
tion, namely K > a — 1, is derived by considering the 
other templates in the chain. In fact, this very simple 
result quantifies nicely the notion that the cooperative 
couplings must reach a certain minimum strength so as 
to balance the competition between templates. 

The analysis of the hypercycle, i.e. to = n, is a little 
more involved. ^From x 2 = we get $ — Q = KQx\ 
which, inserted in the equations x 3 = . . . = x n = 0, 
yields x\ — x 2 = ■ ■ ■ = x n -\ and 



x k 



x k (Q + KQx^x - $) 



(10) 



(11) 



which does not contain x\. Again, in the steady-state 
regime the first equation yields $ = Q, implying KQxi = 



Xi 



K 



(15) 



Finally, using these results in Eq. (||) we find that x\ is 
given by the roots of the quadratic equation 



nKx\ - {KQ 



l)si + l-Q = 0. 



(16) 
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For K < (a — 1) /in, this equation has real roots for all 
Q > 0, otherwise it has real roots for Q > Qh where Qh 
is the unique positive root of the equation 



K 2 Q 2 h + 2K (a - 1 + In) Q h + (a - l) 2 - inK = 0. 



(17) 



In particular, for large K we find Qh ~ 2\/n/K. Further- 
more, it can be easily seen from Eqns. ([l^) and ( |l6| ) that 
x n vanishes at Q = Q n with Q n given in Eq. ( p"4| ) . To un- 
derstand the role of Qh and Q n (we note that Q n > Qh) 
in delimiting the region of existence of the mcoexistence 
state we must look at the behavior of the two real roots 
of Eq. (|l6|). Let us denote them by xf and x^ with 
x i > Xi , which, according to Eq. dl5|), correspond to 
x^ and x~, respectively. Of course, these roots become 
identical at Q = Qh and so the two solutions for x n will 
vanish simultaneously only at the value of K — Kh at 
which Qh equals Q n . Explicitly, we obtain 



K h = {a-l) [n (o+l) - 1] 



(18) 



by inserting Eq. (|l4| ) into Eq. (|l7j ) . Although both roots 
x^ and are in the simplex (0,1), this is not so for 
x+ and x~. In particular, for K < Kh both concen- 
trations are negative within the range Qh < Q < Qn- 
However, while x+ becomes positive for Q > Q n (it 
vanishes at Q n ), x~ remains always negative. Since 
K h > (a - iy /in the same conclusion holds in the 
range K < (a - l) z /in as well, provided we define 
Qh = in this region. The situation is reversed for 
K > Kh- both concentrations are positive within the 
range Qh < Q < Qn, but now it is x~ that vanishes at 
Q n and becomes negative for Q > Q n while x+ remains 
always positive. Despite the small region in the parame- 
ters space where the root x\ yields concentrations inside 
the simplex, the linear stability analysis discussed in the 
sequel indicates that this solution is always unstable, so 
we only need to consider the root x\ . Thus the range of 
existence of the hypercycle fixed point m = n is Q > Q n 
\iK<K h and Q > Q h if K > K h . 

In models without error-tail, i.e., pure replicator equa- 
tions a much stronger statement on coexistence is possi- 
ble. The "time average theorem" |^(| states that if there 
is a trajectory along which a certain subset of templates 
J survives, then there is a fixed point with exactly the 
J-coordinates non-zero. While we have not been able to 
prove the "time average theorem" in full generality for 
Eqns. (Q) and (||), it is easily verified for free chains. 
Hence, if there is no m-coexistence equilibrium, then 
there is no trajectory at all along which the templates 
I\ through I m survive. 

Hitherto we have determined the ranges of the pa- 
rameter Q where the m-coexistence states are physically 
meaningful, in the sense that Xi € (0, 1) Vi. The next 
step is to find the regions where these states are locally 
stable. 



IV. STABILITY ANALYSIS 

In order to perform a standard linear stability analysis 
of the fixed points obtained in the previous section, it is 
convenient to rewrite the kinetic equations (jl]) and (|^) as 
follows 



Xi XiFi (x) % — 1, . . . , n 



(19) 



where 



Fi(x) = AiQ + KQxi^i -A e - J2 x i {Aj-A e +K Xj -i) 

3 

(20) 

and we have used the constraint (|j) to eliminate x e . The 
stability of a fixed point is ensured provided that the 
real parts of all eigenvalues of the n x n Jacobian J are 
negative. In our case the elements of the Jacobian are 
given by 



J ij &ij Ff 



m . 

2-i n J t, 

OXj 



(21) 



The evaluation of the eigenvalues is simple only for the 
zero-coexistence state, since in this case the Jacobian is 
diagonal with elements Jn = aQ — 1 and Ja — Q — l, i > 
1. Therefore this steady state becomes unstable for 
Q > 1/a, which coincides with the lowest replication 
accuracy required for the existence of the 1-coexistcncc 
state. However, for a general m-coexistence state we have 
to resort to a numerical evaluation of the Jacobian eigen- 
values. 

Fortunately, in the case of chains < m < n there is an 
alternative way to look at the stability of the fixed points, 
as hinted by the stability analysis of the zero-coexistence 
state, which becomes unstable due to the emergence of 
the 1-coexistence state. In fact, it can be easily seen 
that any perturbation of the m-coexistence fixed point 
which makes the concentration x m +i non-zero will be 
amplified if A m+ iQ + KQx m — $ is positive. For m > 
we use $ = aQ and A m+ i = 1 together with the value of 
x m given in Eq. (13) to obtain the following (necessary) 
condition for the stability of the m-coexistence state, 



Q < Q m +i m > 0, 



(22) 



with Q m given in Eq. (14). Hence the maximum value 
of Q allowed for the stability of the m-coexistence state 
coincides with the minimum Q required for the existence 
of the (m + l)-coexistence state. Interestingly, though 
for m = we have $ = 1, Ax = a and Xq = 0, condition 
( p2| ) holds true in this case too. 

At this point two caveats are in order. First, the entire 
argument leading to the stability condition (g2|) is flawed 
if the (m + l)-coexistence state happens to be unstable. 
Therefore, we must guarantee via the numerical evalu- 
ation of the Jacobian eigenvalues that the Z-coexistence 
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state is stable before using that condition to study the 
stability of chains with m < I. In particular, we have 
carried out the numerical analysis for the hypercycle so- 
lution I — n and found the following general results: 

(i) For n < 4, it is always stable; 

(ii) for n = 5, it is stable in a very small range of Q 
above Q5; and 

(iii) for n > 6, it is always unstable. 



0.3 



Second, the derivation of the stability condition (22) is 



based on the analysis of a single eigenvalue of the Jaco- 
bian and so it does not yield a sufficient condition for 
the stability of the fixed points. Nevertheless, we have 
verified through the numerical evaluation of all n eigen- 
values that, provided the (m+ l)-coexistence state is sta- 
ble, the eigenvalue associated to fluctuations leading to 
an increase of the chain length is the first one to become 
positive. 
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0.1 




(0) 



V. DISCUSSION 

Combining the existence and the stability results de- 
rived in the previous sections we can draw the phase- 
diagrams in the plane (K, Q) for fixed a and n. In par- 
ticular, for n < 4 the m-coexistence state is stable within 
the interval 



Qm<Q <Q 



m < n 



(23) 



with Q m given in Eq. (14), provided that K > a — 1. 
Interestingly, for fixed Q, Eq. (^3|) shows that the incre- 
ment 5K in the catalytic coupling needed to incorporate 
a new template into the chain is 



SK = 



(a -if 
aQ- 1 ' 



(24) 



regardless the number of elements in the chain. The case 
K < a — 1, for which no chains with m > 1 are allowed, 
does not require any special consideration. In fact, we 
find that the only stable states are the zero-coexistence 
(Q < 1/a) and the 1-coexistence (1/a < Q < 1) states. 
However, since Q2 < 1 only for K > a — 1, this result is 
consistent with Eq. (|23|). 

The n-coexistence state (i.e., the hypercycle solution) 
is stable for Q > Q n if K < Kh and for Q > Qh, oth- 
erwise, where Kh and Qh are given by Eqns. ( |l8|) and 
(fl7|), respectively. We define the error threshold of the 
hypercycle as the value of the replication accuracy Q that 
delimits the region of stability of the n-coexistence state. 
The phase diagram for a — 10 and n = 4 shown in Fig. ^ 
illustrates the major role played by Kh in the hypercyclic 
organization: only for K > Kh the hypercycle becomes 
more stable than a chain of the same length. Another 
important quantity is the value of K, denoted by JsT c , at 
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FIG. 2. Phase-diagram in the space (K, Q) for n = 4 and 
a = 10 showing the regions of stability of the diverse coexis- 
tence states. The numbers between parentheses indicate the 
number of coexisting templates. Regions of bistability ap- 
pear for K > Kh- The thin lines are (from bottom to top) 
Q\,Q2,Q3 and Qi. The thick line is Qh- 

which Qh equals 1/a, the minimal replication accuracy 
of the quasispecies. It is given by 



K c = a (a - 1) 2n - 1 + \fn (n - 1) 



(25) 



Beyond this value the error threshold of the hypercycle 
Qh is smaller than that of the quasispecies. Moreover, as 
mentioned before, for large K, it vanishes like 1/VK. A 
rather frustrating aspect of Kh and K c is that both are of 
order a 2 , indicating then that the productivity of catalyt- 
ically assisted self-replication is much larger than that 
of non-catalyzed self-replication. While this is obviously 
true for biochemical catalysis, it is difficult to argue for 
the existence of such efficient catalysts in prebiotic con- 
ditions. On the other hand, we can take a different, more 
optimistic viewpoint and argue that modern biochemical 
catalysts (enzymes) are so efficient because their precur- 
sors had to satisfy the stringent conditions imposed by 
surpassing Kh- 

In Fig. Q we present the phase-diagram for n — 5. The 
main difference from the previous figure is that the 5- 
coexistence state is stable only within the thin region be- 
tween Q5 and the dashed curve, obtained through the nu- 
merical evaluation of the Jacobian eigenvalues. As these 
curves intersect at some K < Kh, the 5-membered hy- 
percycle is not very interesting, since it has the same 
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FIG. 3. Same as Fig. ^ but for n = 5. There are no stable 
fixed points above the dashed curve. The solid curves are 
(from bottom to top) Qi, Q2, Q3, Q4, and Q5. 



characteristics of a chain of length m = 5. To confirm 
this result we have carried out the numerical integra- 
tion of the kinetic equations using the ninetieth-order 
Runge-Kutta method. The results are shown in Fig. ||, 
which illustrates the time evolution of the concentrations 
Xi {i = 1, . . . , 5) inside and outside the region of stability. 
Although the behavior pattern in the region of instabil- 
ity seems periodic, we have not explored completely the 
space of parameters to discard the existence of chaotic 
behavior. Hence we use the term complex dynamics to 
label this region in Fig. |3[ We note that the phase- 
diagram shown in this figure describes also the regions 
of stability of hypercycles and chains of size n > 5, since 
m-coexistence states with m > 5 are always unstable. 

An interesting limiting case which deserves special at- 
tention is the symmetric hypercycle (a — 1). According 
to the argument put forward in the beginning of Sec. Ill, 
the only fixed points in this case are the zero-coexistence 
state and the hypercycle, i.e., chains are not allowed. 
Moreover, Eq. (|lfi|) yields x\ = X2 = ■ ■ ■ = x n where x\ 
is given by Eq. (jig) with a replaced by 1. The analysis 
of the roots of that quadratic equation and the numeri- 
cal evaluation of the Jacobian eigenvalues yield that the 
symmetric hypercycle is stable for 



K>^(l-Q), (26) 
provided that n < 4. The region of stability observed 
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FIG. 4. Time evolution of the five concentrations of the 
templates composing a hypercycle of size n = 5 for a = 10, 
Q — 1, and (a) K = 37 (inside the region of stability) and (b) 
K = 40 (outside the region of stability). The initial state is 
Xi(0) = 0.2 VI 



in Fig. U for the 5-coexistence state does not appear 
in the symmetric case a — 1, so it must be a conse- 
quence of the asymmetry in the productivity values of 
the non-catalyzed self-replication reaction. We note that, 
differently from the asymmetric case (a > 1), the zero- 
coexistence state is always stable. 

For the sake of completeness, we present some results 
on the elementary hypercycle (Ai = 0, i = 1, . . . , n) cou- 
pled to an error tail (A e — 1) via the imperfect catalyti- 
cally assisted self-replication. Inserting these parameters 
into Eq. (|o|) and setting Fj = Vt yields X\ — . . . = x n 
with x\ given by the larger root of the quadratic equation 



Knx\ - (n + KQ) xx + 1 = 0, 



(27) 



since we have verified that the smaller root is always un- 
stable. As in the symmetric case discussed above, for 
n < 4 the stability condition coincides with the condi- 
tion for real xi, namely, 



^ „ / n n 

Q > 2W . 

^ - V K K 



(28) 



Thus the term in the right hand side of this inequality 
yields the error threshold of the elementary hypercycle. 
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FIG. 5. Time evolution of the five concentrations of the 
templates composing a free chain of size n = 5. The parame- 
ters and initial state are the same as for Fig. 0(b). 



Before concluding this section, we must note that the 
chains of size to considered hitherto arc bonded in a hy- 
percycle of size n > to. We could study free chains of 
size n as well by simply setting xq = in the kinetic 
equations (|l|) and (^). Not surprisingly, the results are 
essentially the same as for bonded chains, with Q m play- 
ing a similar fundamental role in delimiting the regions of 
stability of the shorter chains (to < n). Although a full 
discussion of the stability of the complete chain (to = n) 
is beyond the scope of this paper, we point out that free 
chains with n > 4 are always unstable. Moreover, as il- 
lustrated in Fig. ||, the oscillatory behavior pattern of the 
template concentrations, which ensues a dynamic coexis- 
tence among all templates in the chain, is similar to that 
observed in the hypercycle (compare with Fig. ^) . In this 
sense, the free chains seem as good as the hypercycle to 
attain that kind of coexistence. However, we must em- 
phasize that, for sufficiently large K (i.e. K > Kh), the 
fixed points describing the coexistence of n < 4 templates 
in a hypercycle are much more robust against replication 
errors than their counterparts in the free chains. 



VI. CONCLUSION 



the hypercycle (Qh) becomes smaller than that of the 
quasispecies {Q\ = I /a) for catalytic couplings (K) of 
order of a 2 , where a is the productivity value of the mas- 
ter template. In particular, Qh vanishes like ^JnjK for 
large K. Perhaps, even more important is our finding 
that the asymmetry in the non-catalyzed self-replication 
reaction (a > 1) entails the existence of chains of size 
n < 5. We note that these chains are unstable in the 
symmetric hypercycle as well as in the fully connected 
network [p^| , ^6| . 

Adding to the scenario for the emergence of the hyper- 
cycle described in Sec. ||, which starts with an isolated 
quasispecies, our results indicate that the chains may well 
be the next step in this complex evolutionary process. In 
fact, according to Eq. (|24|) the strengths of the catalytic 
couplings needed to form a chain are of order a, while the 
hypercycle only acquires its desirable stability character- 
istics for strengths of order a 2 (see Eq. (18)). Although 
we realize that an evolutionary step leading from chains 
to hypercycles is still a major one, it is certainly much 
more plausible than a direct transition from quasispecies 
to hypercycle. In any event, we think that the emergence 
of the hypercycle can be explained as a series of plausible 
smooth transitions, without need to postulating the hy- 
percycle as an unique event in prebiotic evolution. In this 
vein, this work represents a modest first step to tackle 
this fundamental problem within a firm basis. 

To conclude we must mention an alternative resolution 
for the information crisis in prebiotic evolution which 
has received some attention recently pi[| , namely, the 
stochastic corrector model [p2| . This model builds on 
ideas of the classical group selection theory for the evo- 
lution of altruism p3| , since it considers implicative tem- 
plates competing inside replicative compartments, whose 
selective values depend on their template composition. 
However, the chemical kinetics equations governing the 
dynamics of the templates inside the compartments dis- 
play a non-physical (non-integer exponents) dependence 
of growth on template concentrations. It seems to us that 
this basic assumption of the stochastic corrector model 
must be relaxed, or at least justified, before it can be con- 
sidered an important alternative to the more traditional 
approach based on the hypercycle and its variants. 
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Our study of the steady states of an asymmetric hy- 
percycle composed of n error-prone self-replicating tem- 
plates indicates that, for n < 4, the error threshold of 
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